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Dynamics of strongly interacting Fermi gases, consisting of a 50-50 mixture of two 
different fermionic species, is investigated. For the equation of state we consider 
a Pade [2/2] approximations, which gives the weak-coupling perturbative formula 
(up to 4th order) in the low density regime, the unitary-limit Monte Carlo result 
in the high density regime, and reproduces the 4-fermion prediction for dimer- 
dimer scattering length in the BEC region. We use a time-dependent LDA to 
derive various properties of the Fermi gas under a harmonic confinement and 
compare them with the data of very recent experiments of ^Li atoms across a 
Feshbach resonance. 
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In this talk the dynamics of strongly interacting dilute Fermi gases (dilute in the sense that 
the range of interatomic potential is small compared with inter-particle spacing) consist- 
ing of a 50-50 mixture of two different states and confined in a harmonic trap Vexti'r) = 
(m/2)(c<j^(a;^ + y^) + uj'^z'^) is investigated in the single equation approach to the time- 
dependent density-functional theory [1,2] 



where Vxc{r,t) = [ "'g^ ]n=n{r,t), the density of the system is n{f,t) =| ^(f, t) p, and the 
velocity field v(f,t) = ;i(^*V* - W**)/(2zmn(f, t)). 

Let us come back to the variational formulation of the Kohn-Sham time-dependent theory 



where [■^ > is a product of two Slater determinants, one for each internal state built up by 
the Kohn-Sham orbitals ipi, and H — T + U is the LDA Hamiltonian. 
Eq.(l) can be derived from Eq.(2) using two approximations 

(i) local transform ip^ ~ (piCxp{ihx/Tn), where ipi and x ^^^^ I'sal functions, 

and 

(ii) < 0|T|0 >^ J{tTF{n) +twinj)n{f,t)d\ 

where |0 > is the product of two Slater determinants built on (^j alone, tTpip) — 
(37r^)^/^(3/i^/(10m))n(f, i)^/^ is the Thomas-Fermi kinetic energy density, and tw{n) — 
{h^ /{%m)){ynf /n is the original von Weizsacker density (OWD). 

We note here that the approximations (i) and < (t)\T\(t) >^ J tTF{'n)n{f, t)d^r lead to the 
hydrodynamic approximation (HA) [3]. But near the surface the Hartree-Fock (HF) type 
densities are proportional to the square of the last occupied state. Therefore, the OWD is 
important in this case and it is expected to determine the asymptotic behavior of the density 
at large distances. It is also expected that the OWD is important in the case of the tight 
radial trapping, A ^ 1. 

Taking into account that in the hmit a — > —0, where a is the scattering length, Eq.(2) 
leads to exact equations describing scaling properties of a ideal gas trapped in a time- 
dependent harmonic potential [4], we do expect that the Kohn-Sham approach may be 
an appropriate to describe dynamics of cold fermionic gas in the regime when superfluid 
quantum hydrodynamics, Eq.(l) is not applicable. 
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For the negative S-wave scattering length between the two fermionic species, a < 0, in 
the low-density regime, kp \ a 1, the ground state energy per particle , e(n), is well 
represented by an expansion in power oi kp \ a \ [5] 

e(n) = 2Ef[^ - 7^A;f | a \ +0.055661(A;f \ a \f - 0.00914(A;f | a |)^ + -], (3) 
10 on 

where Ep = fi^kp/{2m). 

In the opposite regime, a — > — oo ,the Bertsch many-body problem, e{n) is proportional 
to that of the non-interacting Fermi gas 

where a universal parameter /3 [6] is estimated to be /3 = —0.56 [7]. 

In the a —>■ +0 limit the system reduces to the dilute Bose gas of dimers [8] 

e(n) = EFi-l/ikpaf + a„A;f/(67r) + ...), (5) 

where Um is the boson-boson scattering length. Solution of 4-fermion problem for contact 
scattering provided the value ~ 0.6a [9]. 

In Rcfs.[l,2] it has been proposed a simple interpolation of the form e(n) ~ EpPikFO,) 
with a smooth function P{x) mediating between the known limits. 

For the negative a it has been proposed a [2/2] Fade approximant for the function P{x) 

where 5i = 0.106103, 82 = 0.187515, ^3 = 2.29188, 5^ = 1.11616. Eq.(6) is constructed 
to reproduce the first four terms of the expansion (3) in the low-density regime and also 
to reproduce exactly results of the recent Monte Carlo calculations [7], {3 = —0.56, in the 
unitary limit, kpo, — ^ —00. 

For the positive a case ( the interaction is strong enough to form bound molecules with 
energy Eynoi) it has been considered a [2/2] Fade approximant 

^' 2Ef^ l + a^x + a^x-'' ^> 
where parameters a are fixed by two continuity conditions at large x, 1/x — > 0, and by 
two continuity conditions at small x. For example, a\ — 0.0316621, a2 — 0.0111816, — 
0.200149, and = 0.0423545 for = 0.6a. 
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The aspect ratio, presented in Fig. 1, shows that the effect of inclusion of the OWD 
(quantum pressure) on the expansion of superfluid for the conditions of Ref.[6] is about 1%. 

Fig. 2 and Fig. 3 show the comparison between [2/2] Fade approximations, Eqs.(6,7), 
and the lowest order constrained variational (LOCV) approximation [10] and the BCS mean- 
field theory [11] for e{n). The LOCV calculations agree very well with the [2/2] Fade 
approximation results on the BCS side (a < 0). It is evident the difference between our 
results and the BCS mean-field theory calculations. For example, the BCS mean-field gives 
/3 = -0.41. 

The predictions of Eq.(l) with e(n) from Eq.(6) for the axial cloud size of strongly inter- 
acting ^Li atoms are shown in Fig 4. It indicates that the TF approximation of the kinetic 
energy density is a very good approximation for the experimental conditions of Ref.[12], 
A^A ^10^ (inclusion of the OWD gives a negligible effect, < 0.5%). 

In Fig. 5, we present the calculations for the frequency of the radial compression mode 
(jJrad as a function of the dimensional parameter (N^^^a/aho)^^ in the case of an anisotropic 
trap {cUx — cUy — u!±, uJzl^i. — One can easily see that the corrections to the hydrody- 
namic approximation (HA) (inclusion of the OWD) are important even for relatively large 
N and XN . For example, the correction to ojrad in unitary hmit is larger than 11% and 25% 
for A = 10-2, N = 10^ and A = 10"^ N = 10^ respectively. 

In the HA, lo^ad is independent of N for a fixed {N^^^a/ aho)~^ ■ The deviation from 
this behavior does not demonstrate the cross-over to the ID behavior, since AA^ > 1. It 
demonstrates that the validity of the HA depends on the properties of the trap. We note 
here that the collective modes of the Fermi gas under harmonic confinement in the framework 
of the hydrodynamic approximation was considered recently in [13,14]. 

In Fig. 6, the calculated radial compressional frequency is compared with experimental 
data [15] in the BCS-BEC crossover region. There is a very good agreement between calcu- 
lations and experimental data [15]. 

However our calculations for uj^-oa disagree with experimental data of Ref.[16]. It is 
well known that the hydrodynamic equation is expected to be applicable for describing the 
macroscopic excitations of the system up to energies of the order of the energy gap. A, needed 
to break-up a Cooper pair. But for the trapped gas A is a function of r (A decreasing when 
we go away from the center) . It is naturally to assume that condition of the applicability of 
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hydrodynamics to describe the macroscopic excitations of the system at T = is 

^ « 1, (8) 

where fl is the frequency of the macroscopic excitations, A = / n{r)A{r)df/N. To calculate 
A we have used results of Refs.[17,18]. 

Fig. 7 shows the comparison between the average energy gap and frequencies of the trans- 
verse and axial breathing modes. It may explain why our calculations disagree with exper- 
imental data of Ref.[16], see also Ref.[19]. Taking into account the trap difference between 
Ref. [15] and Ref.[16] we expect that the Duke University group will reproduce the 910G 
strong change in the Qrad of Rcf.[16] at B ^ lOOOG. 

As for axial mode, interesting results may start at i? ~ 1250G if ujho of the trap is about 
27r60 Hz. 

We thank J.E. Thomas for his interest and providing us with the experimental data. 
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Fig. 1. Aspect ratio of the cloud of the = 7.5 x 10'^ ^Li atoms as a function of time after 
release from the trap {uj± = 27r x 6605Hz, cj^ = 27r x 230Hz). The circular dots indicate 
experimental data from the Duke University group [6]. The solid line and the dashed line 
represent theoretical calculations in the unitary limit (a — oo) including the quantum 
pressure term and in the hydrodynamic approximation, respectively. 
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Fig. 2. The ground state energy per particle, e(n), in units of 3h^kp/{l0m) as a function 
of the gas parameter {kpa)'^. The sohd hne, the long dashed line and the short dashed 
line represent the results calculated using the [2/2] Pade approximation, Eq.(6), the LOCV 
approximation, and the BCS mean-field theory, respectively. 
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Fig. 3. The ground state energy per particle, e(n) + \Emoi\/'^-i in units of h'^kp/{2m) as 
a function of the gas parameter {kpa)~^. The dashed hne, the dotted-dashed hne and the 
sohd hne represent the resuhs calculated using the BCS mean-field theory, the [2/2] Fade 
approximation, Eq.(7), with = 2a, and am = 0.6a, respectively. 
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Fig. 4. Axial cloud size of strongly interacting ^Li atoms after normalization to a non- 
interacting Fermi gas with A?" = 4 x 10^ atoms as a function of the magnetic field B. The 
trap parameters are uj± = 27r x 640Hz, uj^ = 2it{Q00B /kG + 32)^/^Hz. The solid line and 
dashed line represent the results of theoretical calculation that includes the OWD or uses 
the TF approximation for the kinetic energy density, respectively. The circular dots indicate 
experimental data from the Innsbruck group [12]. 



9 



o 
c 

<D 
0) 



CO 

c 
o 

0) 

<D 
i_ 

Q. 

E 
o 
O 

OS 
DC 



2.4 



2.3 



2.2 



2.1 



1.9 



1.8 



1.7 




■10 



-8 



-6 



-4 



-2 







X 



Fig. 5. Radial compressional frequency, corad, of the cloud of the N — 10^ fermions (solid 
line) and = 10^ fermions (dashed line) in unit of uj± as a function of the dimensional 

parameter 'X.= (N^^^a/aho)^^- The trap parameter A is assmned to be equal to 10~^. The 
lower line (dashed-dotted line) represents the results in the hydrodynamic approximation, 
in which iVrad is independent of N for a fixed {N^^^aj aho)~^ . 
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Fig. 6. The radial compressional frequency as a function of y^={N^^^a/ aho)^^ ■ The sohd 
hne and the dashed hne represent the resuhs calculated using the [2/2] Fade approxima- 
tion with ttm = 0.6a and = 2a, respectively. The solid circles with error bars are the 
experimental results given by the Duke University group [15]. 
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Fig. 7. The average energy gap in units of huJhoN^^^ of an elongated trapped = 4 x 10^ 
Fermi atoms (A = 0.045) as a function of the parameter X = {kp{0)a)~^ (sohd hues). The 
dashed hues in the main plot and in the inset are the frequencies of the transverse and axial 
breathing modes in the same units (hcUhoN^^'^), respectively. 
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